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SUMMARY 


Optical  wavelength  radiation  from  the  wakes  of1  small  gun- launched  models 
(spherically  blunted  bodies)  has  been  measured  at  speeds  from  5  to  10  km/sec 
and  at  stream  densities  from  0.001  to  0.2  atmosphere.  Results  indicate  that 
radiation  in  the  wavelength  range,  0.2  to  1.0  micron,  is  caused  primarily  by 
the  presence  of  ablation  material  in  the  wake.  The  radiation  intensity  and 
spectral  quality  differ  considerably  among  the  various  plastic  materials 
tested,  (a  polycarbonate;,  polyethylene,  an  epoxy,  and  polymerized  formaldehyde). 

For  conditions  at  which  a  majority  of  the  data  were  obtained,  the  peak 
radiant  intensity  from  the  wake  varied  directly  with  stream  density  and 
approximately  with  the  7.5  power  of  the  flight  velocity ,J 

The  variation  of  intensity  of  the  wake  radiance  with  downstream  distance 
along  the  wakes  of  polycarbonate  models,  following  the  region  of  peak  inten¬ 
sity,  could  be  approximated  by  exponential  decay  functions.  The  decay  rates 
varied  with  the  square  root  of  the  stream  density. 

When  the  intensity  of  the  wake  radiation  is  high,  the  radiative  heat 
load  on  the  body  base  can  be  significant.  Results  of  calculations  of  heating 
rates  on  the  model  base  are  presented  and  relations  for  scaling  heating  rates 
are  discussed. 


IRTROBUCTION 


Thermal,  visible,  and  ultraviolet  radiation  from  flow  fields  of  hyper¬ 
velocity  bodies  has  been  extensively  studied  recently  because  of  its  contri¬ 
bution  to  the  heating  of  high-speed  entry  bodies  and  because  of  its  possible 
use  in  tracking  entry  vehicles.  Measurement  of  the  radiative  intensity  and 
its  geometric  distribution  has  also  been  a  valuable  aid  in  the  understanding 
of  various  details  of  the  flow  around  hypersonic  bodies  (ref.  l) .  Because 
some  evidence  indicates  that  the  occurrence  of  bright  luminous  wakes  is 
associated  with  a  dense  electron  population,  information  showing  intensity 
and  distribution  of  wake  radiance  may  be  useful  in  conjunction  with  microwave 
scattering  experiments  (ref.  2). 


Although  some  information  on  the  spectral  and  geometric  characteristics 
of  wake  radiation  from  large-scale  ablating  bodies  can  be  obtained  from 
photographic  and  radiometric  observations  of  their  flight  into  the  atmosphere, 
it  is  not  economically  feasible  to  conduct  large  numbers  of  such  experiments, 
and  when  they  are  made,  results  can  be  difficult  to  evaluate  because  of  the 
necessarily  large  distances  of  the  instrumentation  from  the  entry  trajectory. 

Many  characteristics  of  radiating  wakes  can  not  be  predicted  well  from 
theory,  because  such  predictions  require  the  solution  of  a  combination  of 
complex  problems:  determination  of  transition  locations,  chemistry  and 
reaction  rates  of  high  enthalpy  reactions,  ionization  effects,  and  radiative 
properties  of  the  species  involved.  Some  wake  characteristics  have  been  the 
subject  of  theoretical  investigations  by  various  authors.  References  3,  4, 
and  5  present  studies  aimed  at  the  prediction  of  thermal  properties,  pres¬ 
sures,  degree  of  ionization,  and  the  effect  of  distance  on  the  lateral  extent 
of  the  wake  region,  and  reference  3  includes  results  of  wake -radiation  cal¬ 
culations  for  a  nonablating  body. 

During  a  test  program  initiated  at  Ames  Research  Center  to  study  radia¬ 
tive  heating  from  shock  layers  of  small  gun- launched  models  in  flight  at 
speeds  ranging  from  5  to  10  km/sec  and  at  free-stream  densities  from  0.001 
to  0.2  atmosphere,  the  wake  radiation  was  also  measured.  The  shock-layer 
results  have  been  described  in  references  6  and  7,  which  report  studies  of 
radiation  from  air  in  the  gas  cap,  and  references  8  and  9,  which  present 
results  of  tests  to  determine  the  effects  on  radiative  heating  of  ablation 
products  in  the  gas  cap. 

The  purpose  of  the  present  report  is  to  present  and  correlate  the  wake 
observations.  Predictions  (based  on  the  experimental  data)  of  scale  effects 
on  radiative  heating  of  the  afterbody  and  on  spectral  power  emitted  by  wakes 
are  included. 


SYMBOLS 

Ap  effective  ablation  area,  cm2 

Av  cross-section  area  of  viscous  wake,  cm2 
B  stagnation-point  velocity  gradient  factor 

D  body  diameter,  cm 

Heff  effective  heat  of  ablation,  watt  sec/g 

I  radiant  power  per  unit  volume,  watts/cm3 

lb  integral  appearing  in  equation  (A15) 

K  exponential  decay. factor 


2 


m 

Mb 

K 

p 

p 

Ps 

q. 

R 

% 

R 

T 

u 


W 

wA 

X 

z 

r 

5 

A 

M* 

I 

P 

a 


ablation  mass  rate,  g/sec  cm2 
mass  flow  in  boundary  layer,  g/sec 
mass  flow  in  wake,  g/sec 
pressure,  mm  Hg 
radiative  power,  watts 
similarity  parameter 
radiative  heat  flux,  watts/cm2 
body  radius  of  cross  section,  cm 
nose  radius,  cm 
gas  constant 
temperature,  °K 

local  velocity  parallel  to  surface,  km/sec 

free-stream  velocity,  km/sec 

velocity  normal  to  surface 

velocity  of  wake  relative  to  body,  km/sec 

radiative  power  from  a  length  of  wake  equal  to  the  body  diameter,  watts 
spectral  power,  watts/micron 

axial  distance  of  point  in  wake  behind  model  base,  cm 
compressibility 

factor  defined  by  equation  (A3O) 
exponent  defined  by  equation  (A30) 
wavelength,  microns 

absolute  coefficient  of  viscosity,  g/ cm  sec,  also  microns 
exponent  defined  by  equation  (A30) 
density,  g/cm3 
scale  factor,  Rl/RM 
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Subscripts 


b  base  of  body 

e  edge  of  boundary  layer 

g  forebody  shock  layer  (gas  cap) 

L  large  scale 

M  small  or  model  scale 

P  ablation  product 

w  wake 

W  wall 

o  at  location  of  peak  radiance 

co  free  stream 


EXPERIMENTAL  APPARATUS 


Facility 

The  tests  were  conducted  in  the  Pilot  Hypervelocity  Free-Flight  Facility 
of  Ames  Research  Center.  The  facility,  shown  schematically  in  figure  1,  is 
a  countercurrent  device  in  which  7.1-mm-diameter  models  are  launched  with  a 
light -gas  gun  through  an  instrumented  shock -tube  wind-tunnel  test  section 
containing  either  an  opposing  air  stream  having  a  Mach  number  of  6;  or  still 
air.  Reservoir  conditions  for  the  test -section  nozzle  are  provided  by  a 
cold-helium-driven  shock  tube  operated  under  tailored  interface  conditions. 
Most  of  the  data  presented  herein  were  obtained  under  still-air  conditions  at 
velocities  up  to  8  km/sec.  The  1.8  km/sec  opposing  air  stream  was  used  for 
higher  velocities,  approximately  to  10  km/sec.  References  6  and  10  describe 
the  characteristics  of  the  facility.  Model  velocity,  attitude,  and  position 
were  determined  with  conventional  ballistic -range  instrumentation  -  spark 
shadowgraphs  and  electronic  counters.  Properties  of  the  hypersonic  air  stream 
were  determined  from  measurements  of  shock  velocity  and  transient  pressures  in 
the  shock  tube  and  static  pressures  on  the  test -section  wall. 
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Radiometers 


Radiation  was  measured  with  an  array  of  8  to  10  radiometers  consisting 
of  photomultipliers  and  optical  filter  combinations  with  spectral  responses, 
as  shown  in  figure  2  and  as  given  in  the  table  below. 


Wavelength  of  Half -amplitude 


maximum  sensitivity,  band  pass, 

micron  _  micron 


O.235  0.049 

.265  .04 9 

.312  .056 

.369  -°38 

.460  .026 

.518  .059 

.545  .060 

.640  .077 

.758  .071 

.920  .175 


The  radiometers  were  mounted  around  the  test  section  at  right  angles  to  the 
model  flight  path,  as  illustrated  in  figure  1.  Figure  3  shows  the  geometric 
arrangement  of  a  typical  radiometer,  including  the  slit  pair  which  determines 
the  field  of  view  in  the  direction  of  flight.  The  voltage  outputs  of  the 
photomultipliers  are  displayed  on  oscilloscopes  and  photographed,  providing 
a  record  of  the  radiation  as  a  function  of  time. 


Cameras 

Direct  photographs  of  the  luminous  flow  were  obtained  with  image - 
converter  cameras  having  S-ll  photocathode  responses  and  with  a  camera  employ¬ 
ing  a  Kerr -cell  shutter.  Exposure  times  were  typically  0.05  psec.  The  photo¬ 
graphs  aided  in  the  interpretation  of  some  of  the  radiometric  data  and  were 
useful  for  detemining  the  spatial  extent  of  the  luminosity  and  the  presence  of 
turbulence  of  the  wake. 


Models 

The  models  were  made  of  the  following  plastic  materials:  Lexan  (a 
General  Electric  Company  polycarbonate),  high-density  polyethylene,  G.  E.  124 
(an  epoxy  material  designed  to  form  a  char) ,  Teflon  (a  DuPont  polytetrafluoro- 
ethylene) ,  and  Delrin  (a  polymerized  formaldehyde  also  made  by  DuPont) .  Most 
of  the  models  had  spherical-segment  forward  faces  of  0.5  cm  radius  and  short 
cylindrical  afterbodies,  as  depicted  in  figure  4.  The  Teflon  models  had 
tapered  afterbodies  so  that  they  could  be  launched  in  sabots,  because  without 
a  sabot  Teflon  was  not  strong  enough  to  withstand  the  launching  accelerations. 
Sketches  of  a  Teflon  model  and  a  sabot  are  also  shown  in  figure  4. 
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DATA  REDUCTION 


The  total  radiant  flux  from  the  source  per  unit  wavelength  (watts  per 
micron)  per  unit  length  of  the  wake,  in  the  wavelength  band  to  which  the 
radiometer  is  responsive,  was  deduced  from  the  radiometer  outputs.  The  rela¬ 
tions  between  the  radiometer  outputs  and  the  source  power  per  micron  were 
obtained  by  calibration,  as  described  in  reference  6.  The  axial  length  of 
the  region  viewed  by  each  radiometer  was  computed  from  the  known  width  and 
geometric  arrangement  of  the  viewing  slits,  and  with  this  dimension  the  power 
emitted  by  a  length  of  wake  equal  to  the  body  diameter  was  calculated.  Most 
of  the  radiometer  data  are  presented  in  terms  of  this  length. 


RESULTS 


A  typical  set  of  radiometer  records  is  shown  in  figure  5(a) .  The  signal 
produced  by . the  gas-cap  region  of  the  model,  as  it  moved  into  and  out  of  view, 
is  readily  identified  as  a  short  (2  qsec)  initial  pulse.  It  is  followed  by  a 
signal  of  longer  duration  representing  the  wake  radiation  variation  with  time 
(i.e.,  distance)  behind  the  model,  first  increasing  (in  this  example),  and 
then  more  gradually  decaying.  In  some  tests  the  wake  radiation  was  at  a 
maximum  immediately  behind  the  gas  cap  and  then  decreased  more  or  less 
continuously. 

Figure  5(b)  shows  a  typical  image- converter  photograph  of  the  luminous 
regions  around  a  model  in  flight.  In  this  example,  the  material  is  polycar¬ 
bonate,  velocity  is  6.6  km/sec,  and  density  is  O.O72  times  atmospheric. 


Spectra  of  Wake  Region  of  Maximum  Luminosity 

The  oscilloscope  data  from  the  radiometers  were  analyzed  to  determine  the 
spectra  of  the  wake  radiation  at  the  location  of  maximum  brightness.  When 
the  oscilloscope  trace  did  not  exhibit  a  distinct  peak  in  the  wake,  the  max¬ 
imum  was  assumed  to  be  the  value  recorded  at  the  time  just  after  the  bow 
shock  layer  went  out  of  the  field  of  view. 

El  model  material,-  Examination  of  the  intensity  and  spectral 
characteristics  of  the  peak  wake  radiance  indicates  that  the  model  material 
is  of  major  importance.  This  is  shown  in  figure  6,  which  presents  a  com¬ 
parison  of  spectra  obtained  with  models  made  of  polycarbonate,  polyethylene, 

Cr.  E.  124,  polyformaldehyde,  Teflon,  and  aluminum.  The  observed  spectral 
radiation  into  4jt  steradians  in  watts  per  micron  per  model  diameter  is  plotted 
as  horizontal  lines,  the  lengths  of  which  indicate  the  50-percent  band-pass 
widths  of  the  radiometers;  a  curve  is  then  faired  through  these  lines  to 
approximate  the  spectrum. 


It  is  evident  that  the  wakes  of  the  polycarbonate  models  were  much 
brighter  than  those  for  models  made  of  the  other  materials.  For  three  of  the 


materials  shown,  a  large  fraction  of  the  emitted  power  within  the  response 
range  of  the  radiometers  is  in  the  infrared  part  of  the  spectrum. 

Caution  must  he  exercised  when  comparing  some  of  the  results  since  test 
velocities  were  not  all  the  same.  The  effect  of  correcting  for  the  spread 
in  velocities  would  he  to  exchange  the  relative  positions  of  the  polyformalde- 
hyde  and  Teflon  data,  for  example . 

Figure  6  shows  that  with  aluminum,  Teflon,  and  poly formaldehyde  models, 
the  wake  radiation  was  of  very  low  intensity.  As  a  result,  the  accuracy  of 
the  measurements  is  poor  compared  with  that  when  the  wakes  were  bright . 
Inaccuracies  are  due  to  the  difficulty  of  specifically  identifying  the  wake 
contribution,  as  distinguished  from  the  reflections  from  the  test-section 
walls  of  the  much  brighter  gas-cap  radiation.  Although  the  magnitude  of  these 
additional  radiometer  signals  is  not  known,  the  effect  is  small  for  bright 
wakes  but  can  be  significant  when  the  wake  intensity  is  very  low. 

Figure  6  indicates  that  the  aluminum  model,  the  velocity  of  which  was 
6.2  km/sec,  had  a  brighter  wake  than  the  polyformaldehyde  model.  Calculations 
have  indicated  that  transient  aerodynamic  heating  of  the  aluminum  model  after 
launching  would  not  be  sufficient  to  initiate  ablation  (see,  e.g.,  ref.  ll) . 

A  possible  conclusion  from  this  is  that  (if  in  fact  the  model  were  not 
ablating)  the  observed  wake  is  the  result  of  air  radiation.  If  this  is  true, 
it  implies  that  the  ablation  of  the  polyformaldehyde  models  attenuates  the 
air  radiation.  In  general,  however,  these  results,  from  tests  in  which  the 
wake  radiance  was  very  weak,  must  be  considered  as  tentative  because  of  the 
poor  signal- to- noise  ratio  of  the  measurements . 

The  data  in  figure  6  for  the  three  materials  producing  low  wake  intensi¬ 
ties  are  probably  only  useful  to  indicate  upper  limits  of  spectral  power  at 
a  given  wavelength  since  sources  other  than  the  wake  may  be  significant. 

Effect  of  stream  density  on  peak  spectra.-  The  effect  of  varying  the 
free-stream  density  from  0.001  to  0.l8  atmosphere  on  the  spectra  of  the  peak 
wake  radiation  observed  behind  polycarbonate  and  polyethylene  models  is  shown 
in  figure  7-  When  the  spectra  are  compared,  the  fact  that  the  model  veloci¬ 
ties  differ  (ranging  from  5-9  to  7-^  km/sec)  must  be  kept  in  mind,  as  men¬ 
tioned  before.  At  the  higher  densities,  a  major  part  of  the  radiant  power  is 
in  the  infrared  but  with  decreasing  density  it  shifts  toward  the  visible.  At 
the  lowest  density,  0.001  atmosphere,  for  which  data  were  obtained  only  with 
polyethylene,  the  fraction  in  the  infrared  again  increases,  resulting  in  a 
relatively  uniform  distribution  throughout  the  wavelength  range . 

The  previously  mentioned  uncertainty  (in  the  data  for  wakes  that  show 
low  radiation)  resulting  from  the  effects  of  bright  gas  caps  influencing  the 
radiometer  readings,  also  is  present  in  the  data  of  figure  7  &t  short  wave¬ 
lengths.  At  these  wavelengths,  the  wakes  have  relatively  low  radiant  power, 
whereas,  the  gas-cap  power  is  high  (ref.  6).  Thus,  the  uncertainty  is  a 
larger  fraction  of  the  measurements  in  the  ultraviolet  and  short  wavelengths 
than  in  the  visible  and  infrared. 
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Integrated  Peak  Radiative  Intensity 


Effect  of  ambient  air  density,-  The  peak  radiation  spectra  have  been 
integrated  to  obtain  the  total  power  (watts/diameter)  in  the  wavelength  range 
of  the  radiometric  measurements  (0.2  to  1.1  micron).  This  integrated  power 
is  plotted  in  figure  8  as  a  function  of  ambient  density  for  several  veloci¬ 
ties.  Since  the  model  velocities  in  general  differed  from  the  nominal  values 
shown  in  the  figure,  the  data  were  adjusted  using  a  7*5  power  variation  of 
the  radiation  with  velocity.  The  basis  for  choosing  this  power  will  be 
remarked  upon  later. 

Figure  8(a)  presents  the  radiated  power  as  a  function  of  density  for 
three  materials  -  polycarbonate,  polyethylene,  and  G.  E.  124  -  at  a  model 
velocity  of  6*7  km/sec.  The  data  for  all  of  these  materials  show  that  the 
radiant  power  is  essentially  proportional  to  the  ambient  density.  There  is  a 
fairly  large  random  spread  in  the  data  which  is  believed  to  result  primarily 
from  an  inherent  unsteadiness  of  the  wake  luminosity  rather  than  from  a  lack 
of  precision  of  the  measurements.  This  statement  is  based  on  the  fact  that 
the  repeatability  of  the  measurements  of  gas- cap  radiative  intensities  is  much 
better  than  that  of  peak  wake  radiative  intensities.  The  one  data  point  in 
figure  8(a)  for  a  polyf ormaldehyde  model  is  interpreted  as  representing  an 
upper  limit  deduced  from  the  data  as  was  mentioned  before  (cf.  fig.  6). 

Results  obtained  at  higher  velocities,  figure  8(b),  are  generally  similar 
to  those  at  6.7  km/sec. 

Effect  of  velocity.-  In  figures  9(a)  and  9("b)  the  peak  wake  integrated 
radiative  intensity  is  presented  as  a  function  of  model  velocity  in  logarith¬ 
mic  coordinates  for  polycarbonate  and  polyethylene  models  at  several  ambient 
densities.  Included  on  the  figure  is. a  line  representing  a  7.5  power  varia¬ 
tion  with  velocity.  The  lack  of  more  extensive  experimental  data  prevents 
precise  determination  of  velocity  dependence;  however,  at  velocities  up  to 
7  km/sec  for  polycarbonate  models,  and  for  all  the  data  obtained  with  poly¬ 
ethylene  models  (except  perhaps  for  the  lowest  density),  the  wake  intensity 
varies  as  velocity  raised  to  a  power  between  7  and  8.1  At  velocities  above 
7  km/sec,  one  data  point  for  polycarbonate  indicates  a  more  gradual  variation. 


Radiance  Downstream  From  Peak  Region 

Spectra  at  various  distances  behind  the  model,  up  to  3O  model  diameters, 
were  obtained  from  the  oscilloscope  records.  Such  spectra  for  polycarbonate, 
polyethylene,  and  G.  E.  124  models  are  shown  in  figure  10.  Results  shown  for 
a  stream  density  of  O.O72  atmosphere  indicate  little  change  in  spectral 
quality  with  distance;  the  intensities  at  all  wavelengths  decrease  smoothly 
at  about  the  same  rate.  At  a  lower  density,  0.006  atmosphere,  with  polycar¬ 
bonate,  there  was  a  generally  uniform  decrease  with  distance  at  the  longer 
wavelengths  (above  0.55  h)  but  at  short  wavelengths  the  decay  rate  was 

xAn  exponent  of  7.5  was  chosen  for  adjusting  the  data  presented  in 
figure  8  (to  account  for  velocity  variations),  but  choosing  a  different  value 
between  7  and  8  would  not  change  the  results  significantly. 


erratic.  With  polyethylene  and  G.  E.  124  at  low  stream  densities  ( 0.009 
atmosphere  and  0.018  atmosphere,  respectively)  variations  in  the  spectral 
distribution  with  distance  can  be  observed.  With  polyethylene,  this  variation 
consists  of  an  increase  in  the  relative  prominence  of  a  band  near  0.52  \x  and  a 
decrease  in  the  contribution  of  a  bright  band  near  O.37  [x.  With  G.  E.  124, 
the  radiation  at  wavelengths  above  0.5  V-  has  a  considerably  decreased  contri¬ 
bution  at  the  larger  distances. 

The  spectra  obtained  in  the  tests  do  not  have  sufficiently  high  wave¬ 
length  resolution  for  precise  identification  of  the  excited  species  producing 
the  observed  radiation.  However,  figure  10  appears  to  show  some  strong  band 
structure.  Undoubtedly,  an  important  part  of  the  radiation  from  wakes  of 
these  ablating  plastic  models  is  related  to  the  presence  of  a  substantial 
amount  of  carbon  in  some  form  resulting  from  the  thermal  decomposition  of  the 
plastics.  The  spectra  for  polycarbonate  show  a  strong  band  in  the  neighbor¬ 
hood  of  0.5  p,  which  is  probably  attributable  to  the  C2  Swan  system.  The 
intensity  of  this  system  at  low  density  and  a  velocity  of  7.5  km/sec 
(fig.  10(a))  shows  only  a  slight  decay  with  distance  up  to  10  diameters  and 
then  a  much  more  rapid  decay  rate.  At  higher  density,  the  decay  rate  is 
relatively  uniform  up  to  a  distance  of  20  diameters.  For  polyethylene, 
figure  10(b)  indicates  that  this  band  system  is  strong  and  increases  with 
distance  up  to  20  diameters  at  low  ambient  densities,  but  at  a  higher  density 
(O.O72  atmosphere)  it  is  not  prominent  at  any  location. 

An  intense  system  is  evident  in  the  spectra  for  polyethylene  in  the  wave¬ 
length  band  near  O.37  micron.  This  may  be  the  CN  violet  system  (0.388  micron). 
Comparing  figure  10(b)  with  figure  7(b)  shows  that  this  band  intensity  was 
considerably  influenced  by  a  variation  in  velocity,  becoming  quite  weak  with 
a  decrease  in  velocity  from  7*2  to  6.1  km/sec  at  a  stream  density  of  O.O72 
atmosphere.  The  CN  violet  also  appears  to  contribute  importantly  (but  with 
less  relative  prominence  than  for  polyethylene)  to  the  total  power  in  the 
spectra  for  G.  E.  124  (fig.  10(c))  and  at  an  ambient  density  of  0.006  atmos¬ 
phere  to  the  power  in  the  spectra  for  polycarbonate  (fig.  10(a)). 

Throughout  the  range  of  densities  of  the  tests,  the  decay  to  half  ampli¬ 
tude  of  the  luminosity  of  wakes  from  polycarbonate  and  G.  E.  124  models  could 
generally  be  approximated  by  an  exponential  decay  function.  (This  has  also 
been  pointed  out  in  ref.  1.)  A  damping  coefficient  K  has  been  calculated, 
assuming  that  the  decay  of  luminosity  can  be  represented  by  the  relation 

■4  ■  t)  (1) 

W  =  Wie  X  y  (1) 

where  x^  is  the  distance  to  the  start  of  the  decay,  W^_  is  the  radiant  power 
per  unit  length  of  the  wake  at  the  start  of  the  decay,  and  D  is  the  model 
diameter.  From  a  comparison  of  the  damping  coefficients  obtained  from  the 
outputs  of  the  various  radiometers,  it  can  be  determined  whether  the  decay 
rates  depend  on  wavelength.  Figure  11  shows,  as  functions  of  stream  density, 
values  of  K  obtained  from  oscilloscope  records  from  three  radiometers 
responsive  to  the  following,  wavelength  intervals :  a  broad  band  in  the  ultra¬ 
violet  and  visible  (0.21  to  0.6  micron),  a  narrower  band  in  the  visible  (0.48 
to  O.56  micron),  and  a  band  in  the  near  infrared  (0.8  to  1.1  microns).  The 
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three  plotted  curves,  which  each  represent  averages  of  results  from  24  tests, 
indicate  that  at  all  three  wavelength  bands  the  decay  rate  increased  approx¬ 
imately  with  the  square  root  of  the  stream  density  and  that  at  the  different 
wavelengths  the  decay  rates  were  not  greatly  different. 

When  the  radiometer  data  obtained  in  tests  of  polyethylene  models  were 
used  with  equation  (l)  to  calculate  damping  coefficients,  it  was  found  that 
these  coefficients  varied  erratically  both  with  density  and  with  wavelength. 
This  behavior  of  the  decay  rate  is  probably  associated  with  other  anomalies 
that  were  observed  in  the  streamwise  variation  of  the  wake  radiance  of  the 
polyethylene  models,  which  are  discussed  in  the  following  section. 


Variations  in  the  Location  of  Peak  Radiance 

In  tests  with  polyethylene  models,  at  the  highest  and  lowest  densities 
studied,  the  wake  intensity  was  observed  to  be  a  maximum  immediately  behind 
the  location  of  the  gas  cap  and  to  decay  monotonically  thereafter.  At  densi¬ 
ties  between  these  extremes,  there  were  large  variations  in  the  location  of 
peak  luminosity.  These  are  illustrated  in  figure  12,  where  the  wake  radiation 
measured  by  one  radiometer,  sensitive  to  wavelengths  from  0.45  to  0.59  micron, 
is  plotted  as  a  function  of  distance.  The  peak  wake  intensity  is  considerably 
farther  to  the  rear  at  densities  of  0.01  and  0.02  than  at  higher  and  lower 
densities . 

It  can  be  reasoned  that  the  varying  location  of  peak  brightness  may  be 
related  to  a  variation  in  the  characteristics  of  the  viscous  flow  (separated 
boundary  layer)  and  involves  the  rate  of  turbulence  growth.  From  this  reason¬ 
ing,  it  was  believed  that  a  correlation  of  the  data,  based  on  Reynolds  number, 
might  be  possible.  For  this  correlation  to  be  shown  by  experimental  results, 
it  is  necessary  to  include  the  variation  of  more  than  the  single  variable 
density.  Results  from  tests  of  larger  models  (diameter  =  1.52  cm)  were 
available.  These  results  are  shown  in  figure  13  together  with  the  results 
from  the  present  tests,  as  functions  of  Reynolds  number  based  on  free-stream 
conditions.  The  figure  indicates  that  as  the  Reynolds  number  decreased  to 
about  25^000,  the  luminous  peak  moved  rearward.  At  lower  Reynolds  numbers, 
the  peak  is  very  near  the  body.  It  should  be  remarked,  however,  that  for 
these  conditions  it  is  quite  possible  that  a  second  peak  intensity  point 
existed  in  the  wake  at  large  values  of  x/D,  where  the  wake  was  not  scanned. 

Also  included  in  the  figure  are  results  from  reference  12  showing  as  a 
function  of  Reynolds  number  the  location  in  the  wake  at  which  transition  was 
judged  to  occur  on  the  basis  of  the  appearance  of  eddies  in  schlieren  photo¬ 
graphs.  The  similarity  of  conditions  at  which  these  two  phenomena  occur 
suggests  that  they  are  related,  but  further  study  is  required  to  ascertain 
whether  this  apparent  relationship  is  real  or  fortuitous.  The  results 
reported  in  reference  12  were  for  a  much  lower  velocity  -  2.6  km/sec  -  as 
compared  to  6.7  for  the  tests  reported  herein. 

An  image- converter  photograph  obtained  during  the  present  tests  (fig.  14) 
shows  some  of  the  details  of  the  luminous  wake  at  a  Reynolds  number  of  68,400 
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(ambient  density,  0.02  atmosphere).  The  region  in  the  wake  where  the  diameter 
appears  to  increase  as  a  result  of  large-scale  turbulence  coincides  with  the 
location  of  peak  intensity  measured  with  the  radiometer  previously  shown  on 
figure  13. 


Wake  and  Gas-Cap  Radiation  Compared 

The  optical  radiation  emanating  from  the  entire  heated  flow  field  sur¬ 
rounding  a  hypersonic  body  entering  the  atmosphere  is  important  in  locating, 
identifying,  and  tracking  the  body.  In  this  application,  the  photometric 
observations  would  usually  be  made  from  a  great  distance,  so  that,  because  of 
limitations  in  the  spatial  resolving  power  of  the  photometric  instrumentation, 
it  may  be  impossible  to  examine  separately  the  characteristics  of  the  gas  cap 
and  the  wake  regions.  For  this  reason  it  is  useful  to  know  how  the  radiant 
power  from  the  two  regions  compare.  Several  variables  are  involved  in  such 
a  comparison,  including  the  effect  of  the  differences  in  shape  and  extent  of 
the  luminous  regions,  and  the  effect  of  ablation-product  radiation  originating 
in  the  gas  cap.  As  a  preliminary  step,  comparing  the  radiant  power  from  the 
wake  with  the  calculated  power  radiated  from  equilibrium  air  in  the  gas  cap 
gives  a  concept  of  the  relative  contribution  of  the  wake.  The  air  radiation 
values  were  obtained  from  reference  13. 

The  power  from  the  gas  cap  can  be  computed  by  integrating  the  calculated 
local  distribution  of  power  throughout  the  shock  layer,  as  discussed  in  refer¬ 
ence  6,  but  in  the  calculations  discussed  here,  it  was  obtained  by  computing 
stagnation- point  radiation  and  applying  a  correction  factor  to  account  for 
radiation  distribution  in  the  shock  layer.  Values  of  the  correction  factors 
were  estimated  by  the  method  discussed  in  reference  14. 

The  wake  radiation  was  calculated  from  the  integral 

/CO 

E(x) dx  (2) 

where  E(x)  is  power  per  centimeter  of  wake  length.  The  value  of  E  is 
assumed  a  constant  (E0)  for  a  short  distance  (x^)  behind  the  body,  decaying 
exponentially  thereafter.  Then  equation  (2)  becomes 


?w 


E0xi 


-KT(x-Xl) 

e  dx 


(3) 


where  K1  is  the  damping  coefficient,  K,  mentioned  earlier,  divided  by  the 
body  diameter. 

The  ratio  of  the  integrated  wake  output  to  the  output  from  equilibrium 
air  in  the  gas- cap  volume  is  presented  in  figure  15  as  a  function  of  ambient 
density  for  polycarbonate  models  at  two  velocities,  6.7  and  9*7  km/sec.  It 
is  seen  that  the  wake  power  is  much  greater  than  that  of  the  gas  cap,  its 
relative  contribution  becoming  larger  with  decreasing  density. 


11 


Two  factors  should  be  taken  into  account  in  the  interpretation  of 
figure  15;  both  result  from  the  assumption  that  all  the  gas- cap  radiation  is 
the  radiation  from  air  in  thermodynamic  equilibrium.  First,  at  the  lower 
densities  at  which  the  tests  were  conducted,  the  shock-layer  air  was  not  in 
equilibrium  and  measured  shock- layer  radiative  intensities  are  greater  than 
those  calculated  for  equilibrium  conditions.  These  measured  values  are  dis¬ 
cussed  in  reference  6.  The  second  factor  is  that  for  conditions  for  which  the 
wake  intensity  is  high,  there  is  also  a  significant  increase  in  gas-cap  radia¬ 
tion  due  to  the  presence  of  ablation  products  in  the  shock-layer  region.  This 
effect,  which  is  discussed  in  detail  in  reference  8,  would  reduce  the  ratios 
shown  in  figure  15  ut  all  of  the  test  conditions  indicated.  It  should  also  be 
noted  again  that  the  wakes  of  polycarbonate  models  were  the  most  luminous  of 
any  of  those  measured  experimentally. 


APPLICATION  TO  FULL  SCALE 


The  data  and  correlations  presented  thus  far  apply  to  the  small-scale 
models  that  were  tested*  Applying  the  results  to  problems  of  predicting 
radiation  intensity  from  wakes  of  large  bodies  entering  the  atmosphere 
requires  consideration  of  the  effect  of  scale.  Appendix  A  is  an  attempt  to 
predict  peak  intensity  radiation  at  large  scale  on  the  basis  of  present,  data. 
The  following  assumptions  are  made: 

1.  The  radiant  intensity  is  assumed  to  be  a  function  of  only  tempera¬ 
ture,  mass  density  of  the  material  ablated  from  the  body,  and  mass  density  of 
the  air  comprising  the  wake  flow. 

2*  It  is  assumed  that  the  region  from  which  radiation  emanates  is 
uniform  in  temperature,  pressure,  and  chemical  composition  across  each  section 
because  its  components  are  completely  mixed;  these  properties  may  vary  along 
the  wake,  however. 

3.  The  cross-section  area  of  the  viscous  wake  is  assumed  to  be  deter¬ 
mined  by  the  wake  velocity  and  by  the  density  of  the  mixture  of  ablation- 
product  gas  and  air  from  the  boundary  layer. 

A  scaling  relation  is  first  derived  for  bodies  of  different  size  at  a 
given,  constant  free- stream  density.  In  a  second  derivation,  Reynolds  number 
is  held  constant  by  changing  density  in  inverse  proportion  to  scale.  When 
density  is  held  constant  and  the  size  changes,  the  Reynolds  number  varies, 
and  so  the  flow  geometry  is  expected  to  vary.  Which  one  of  these  two  rela¬ 
tions  is  chosen  in  a  particular  calculation  depends  on  whether  the  fluid 
mechanics  or  the  chemistry  of  the  flow  is  expected  to  dominate  the  scale 
effect.  If  the  location  of  boundary  layer  or  wake  transition  and  the  degree 
of  turbulence  in  the  wake  have  significant  effects  on  radiation,  or  if  the 
wake  emissivity  is  so  high  that  the  wake  can  not  be  considered  optically  thin, 
data  obtained  in  tests  conducted  at  full-scale  Reynolds  number  provide  a 
better  basis  for  scaling.  Such  tests  can  also  duplicate  nonequilibrium  flow 
characteristics  in  the  shock  layer,  where  binary  collision  processes  tend  to 
predominate.  Scaling  at  constant  density  permits  the  duplication  of 
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high-enthalpy  chemical  reactions  involved  in  the  boundary  layer  and  wake  flows 
but  does  not  simulate  nonequilibrium  effects  in  either  the  shock  layer  or 
expansion  regions • 


Constant  Density,  Variable  Size 

The  analysis  in  appendix  A  for  the  case  where  free-stream  density  is 
held  constant  indicates  that  the  peak  wake  radiant  power  per  unit  volume 
(watts/cm3)  is  invariant  with  change  in  body  size,  I  =  Im,  while  the  wake 
radius  at  the  location  of  this  peak  radiance  relative  to  the  body  radius 
varies  as  follows,  from  equation  (A25) 

Rw/R  oc  R-1/4 

In  obtaining  these  results  it  was  assumed  that  the  local  pressure  in  the 
wake  was  not  affected.  The  scaling  relation  is  restricted  to  variations  in 
scale  within  which  the  wake  flow  does  not  undergo  a  major  change  in  its  char¬ 
acteristics,  such  as  a  change  from  a  long  laminar  run  to  one  with  laminar-to- 
turbulent  transition  near  or  on  the  body.  An  examination  of  photographs  of 
luminous  wakes  of  models  indicates  that  there  is  a  large  range  of  Reynolds 
numbers  in  which  such  a  major  change  does  not  take  place. 

Experimental  evidence  that  the  relative  wake  diameter  varies  with  the 
-l/4  power  of  the  Reynolds  number  is  shown  in  figure  l6  with  luminous  wake 
measurements  taken  from  photographs  obtained  with  image- converter  cameras  and 
a  Kerr-cell  shuttered  camera. 


Constant  Reynolds  Number 

The  second  analysis  presented  in  appendix  A  furnishes  a  means  of  calcu¬ 
lating  peak  radiative  intensity  of  a  wake  behind  a  large  body  on  the  basis  of 
data  from  tests  of  small  models  at  the  same  Reynolds  number.  The  following 
relation,  expressing  the  effect  of  scale  on  the  peak  radiative  intensity 
(power  per  unit  volume),  is  derived  for  a  flight  velocity  of  6.7  km/sec: 

I,  watts/cm3  oc  R"3/2  (4) 

The  exponent  on  R  was  obtained  from  the  experimental  data;  hence,  it  should 
not  be  expected  to  apply  generally  to  other  materials,  flight  speeds,  etc. 

The  appendix  also  indicates  that  there  is  flow  similarity,  except  for  density- 
dependent  chemical  effects.  The  effect  of  scale  on  the  longitudinal  distri¬ 
bution  of  the  radiation  is  not  known. 
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BASE  HEATING  FROM  WAKE  RADIATION 


The  wake  is  a  source  of  radiant  heating  of  afterbody  surfaces.  In  the 
following  paragraphs,  relations  are  developed  for  this  heating  rate  both  for 
the  models  of  the  present  tests  and  for  full-scale  vehicles. 

The  total  radiant  power  from  an  optically  thin  wake  segment  of  length 
dx  is 

dP  =  itRyl(x)  dx  (5) 

where  l(x)  is  the  power  radiated  from  a  unit  volume  in  the  wake.  The  wake 
is  assumed  to  be  a  line  source  of  radiation,  the  strength  of  which  varies 
along  its  length.  The  average  radiant  flux  intercepted  by  a  body  base  of 
radius  Rb  at  the  distance  x  is 


dt>  =  watts/cm2  (6) 

where  Cl  is  the  solid  angle  intercepted  by  the  base,  (To  obtain  the  actual 
heating  rates,  equation  (6),  of  course,  would  necessarily  include  the  surface 
absorptivity.)  The  solid  angle  is  given  by 


Cl  =  2* 


1  - 


To  calculate  the  heating  produced  by  all  of  the  wake,  equation  (6)  is 
integrated  over  the  wake  length: 


%  = 


r; 


(7) 


Base  Heating  of  Model 

When  equation  (7)  is  used  to  calculate  base  heating,  it  can  be  shown 
that  this  heating  depends  primarily  on  radiation  emitted  from  the  wake  close 
to  the  body  and  that,  among  the  wakes  observed  in  the  tests,  differences  in 
the  longitudinal  distributions  of  wake  radiation  are  of  only  secondary  impor¬ 
tance.  Large  variations  in  this  distribution,  representing  the  entire  range 
of  distributions  observed  in  the  tests  for  which  wake  intensities  were  high, 
had  only  a  small  effect  on  the  heating.  For  this  reason,  a  single  distribu¬ 
tion  function,  which  approximated  the  measured  distributions,  was  used  to 
calculate  base  heating  of  the  model  at  all  stream  densities.  This 
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distribution  consisted  of  a  region  in  which  the  radiation  (watts  per  cm 
length)  increased  linearly  from  a  value  of  0.9  times  the  maximum  at  the  body 
base  to  a  maximum  at  1.3  diameters  behind  the  body,  followed  by  a  region  of 
linear  decrease  to  50  percent  of  the  maximum  at  a  distance  of  10  diameters. 
With  this  distribution  and  the  peak  intensities  shown  in  figure  8,  equa¬ 
tion  (7)  gives  the  base  heating  rates  for  polycarbonate  models  of  0*71  cm 
diameter  tabulated  below.  The  ordinate  W  (fig.  8)  is  related  to  I0,  the 
maximum  value  of  the  parameter  I  (eq.  (7))j  by  the  formula 


Io 


W 


P>0 

W  (fig.  8) 

q-^,  watts/i 

0.003 

1-5 

1.4 

.009 

4.2 

.018 

9.4 

9.0 

.072 

37.2 

35.5 

.180 

90.7 

87.0 

Base  Heating  at  Large  Scale 


In  order  to  examine  the  effect  of  scale  on  base  heating,  equation  (7)  is 
written  in  dimensionless  form  as  follows,  where  R  is  a  characteristic  body 
radius , 


lb 


1  p  P00  (VR)2  ^/R) 

1.0 

x/R 

2  J0  (Rb/R)2 

[(Rb/R)2  +  (x/r)2]1/2 

the  radiant  intensity  as 

/*v 

= 

'l(x/R)1 

w 

xo 

Io  J 

(8) 


and  assuming  Rb  =  R,  equation  (8)  is 


It  has  previously  been  indicated  that  for  the  calculation  of  base  heat¬ 
ing  it  is  necessary  to  know  the  longitudinal  distribution  of  wake  radiation 
in  the  near  wake  only.  The  effect  of  scale  on  this  distribution  has  not  been 
determined;  however,  in  tests  of  the  model,  it  was  observed  that  a  single 
distribution  would  adequately  represent  the  distribution  throughout  a  large 
Reynolds  number  range.  In  the  absence  of  better  information,  and  because  it 
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can  be  shown  that  small  variations  in  the  distribution  do  not  greatly  affect 
the  base  heating  results,  it  is  assumed  that  this  distribution,  described  in 
the  preceding  section,  also  applies  at  large  scale.  Then  the  integral  of 
equation  (9)  is  independent  of  size,  and  the  entire  effect  of  scale  is 
included  in  the  geometry  factors  and  the  value  of  the  peak  radiation  IQ, 
which  can  be  calculated  by  means  of  the  scaling  rules  already  discussed. 

If  the  radiative  heating  of  the  large-scale  body  is  to  be  determined  by 
extrapolating  from  tests  at  the  same  density  and  velocity,  I0  is  constant 
and,  according  to  equation  (A25),  the  wake  radius  varies  with  the  3/4  power 
of  the  scale.  Inserting  these  relations  in  equation  (9)  and  denoting  the 
integral  in  equation  (9)  by  A  gives 

(R3/4\2 

%  °= - 5 -  IoA  K  R1/2IcA  (10) 


ibL  =  °1/24bM  (11) 

If  an  extrapolation  of  data  is  based  on  tests  at  the  same  Reynolds 
number,  at  a  velocity  of  6.7  km/sec,  the  wake  intensity  dependence  is 
expressed  by  equation  (4)  for  the  materials  studied  in  the  present  tests. 

Hence , 


Ip  =  IMa-3' 2 

From  this,  equation  (9)  gives  the  scaling  rule 


%■ 


-  a~1/Z% 


M 


(12) 


Example  of  Scale  Effect 

The  result  of  applying  the  above  scaling  relation  is  illustrated  with 
the  following  example.  Assume  a  scale  factor  of  60,  which  for  constant 
Reynolds  number  corresponds  to  a  body  radius  of  21.3  cm.  For  flight  of  this 
large  body  at  a  stream  density  of  0.003  atmosphere,  the  previously  given 
radiative  heating  rate  that  was  calculated  for  the  model  at  a  stream  density 
of  0.18,  87  watts/cm2,  is  used.  Then  the  radiative  heating  of  the  base  of  a 
large  body  similar  in  shape  to  the  model  is  obtained  from  equation  (12). 


ibL  =  %Ms-1/a  =  (87)  (6o)-1/2 

=  11.2  watts/cm2 
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CONCLUSIONS 


Results  of  ballistic  tests  of  small-scale  ablating  models  at  speeds  up 
to  10  km/sec  have  been  analyzed  and  indicate  the  following  conclusions  rele¬ 
vant  to  the  characteristics  of  their  luminous  wakes. 

1.  Essentially  all  of  the  measured  optical  radiation  from  the  wake  is 
produced  by  the  effects  of  material  ablated  from  the  model.  The  wakes  of 
models  made  of  polycarbonate  were  considerably  brighter  than  those  of  other 
materials . 

2.  The  radiative  intensity  of  the  region  in  the  wake  where  the  bright¬ 
ness  was  a  maximum  generally  varied  directly  with  free-stream  density  at  fixed 
model  velocities. 

3.  The  rate  at  which  the  luminosity  decayed  with  distance  behind  the 
model  varied  approximately  with  the  square  root  of  the  stream  density. 

4.  The  integrated  power  radiated  from  all  the  polycarbonate  model  wake's 
was  about  one  order  of  magnitude  greater  than  that  calculated  for  all  the 
equilibrium  air  in  the  shock  layer  at  the  highest  test  free-stream  densities 
and  about  three  orders  greater  at  the  lowest  test  densities . 

5.  It  was  deduced  from  the  present  results  that  afterbody  heating  from 
wake  radiation  is  mostly  contributed  by  the  near  wake,  usually  that  portion 
of  the  wake  within  ten  body  diameters.  At  constant  stream  density,  the 
heating  rates  are  predicted  to  increase  as  the  square  root  of  the  scale,  and 
at  constant  Reynolds  number  (for  a  velocity  of  6.7  km/sec),  the  heating  rates 
are  predicted  to  vary  inversely  with  the  square  root  of  the  scale. 
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APPENDIX  A 


SCALE  EFFECT  ON  RADIATIVE  INTENSITY 


In  this  appendix,  the  effect  of  scale  on  peak  wake  radiative  intensity  is 
estimated,  on  the  basis  of  calculated  effects  of  scale  on  wake  properties  and 
on  concentration  of  gaseous  ablation  products.  It  is  assumed  that  the  radia¬ 
tion  is  emitted  by  an  inner  viscous  core  containing  the  air  and  ablated  mate¬ 
rial  from  the  boundary  layer.  Effects  of  scale  on  wake  properties  can,  to  a 
large  degree,  be  traced  to  and  correlated  with  effects  in  the  boundary  layer. 
For  this  reason,  the  analysis  will  first  consider  the  boundary- layer  flow. 

Because  most  of  the  ablation  is  from  the  forebody,  the  boundary  layer  on 
only  this  part  of  the  body  is  assumed  to  contribute  significantly  to  the  wake 
core.  Other  assumptions  are  that  (l)  the  boundary  layer  is  laminar;  (2)  the 
wall  temperature  is  constant,  fixed  by  the  sublimation  temperature  of  the  sur¬ 
face  material;  and  (3)  the  ablation  mass  injected  is  gaseous  with  state  and 
transport  properties  identical  to  equilibrium  air.  The  boundary- layer  analy¬ 
sis  is  based  on  solutions  presented  in  references  9  an(^  15*  Reference  15 
shows  that  the  differential  equations  (motion  and  energy)  for  a  two- 
dimensional  stagnation -flow  boundary  layer  with  flow  through  a  porous  wall 
can  be  transformed  so  as  to  yield  similar  solutions  in  the  transformed  system, 
if  the  stream  function  fy  defining  the  flux  through  the  wall  does  not  vary. 
The  transformed  normal  coordinate  is 


=  y(pwuooAwx)1/2 


and  the  transformed  stream  function  is 


1/  2 

fw  =  PwVCMyXUjfy) 


and  up  =  S\|r/3y.  The  coordinates  x  and  y  are  parallel  and  normal  to  the 
wall,  respectively.  An  analogous  solution  for  axisymmetric  flow  is  presented 
in  reference  9*  Making  use  of  the  following  (Mangier  and  Howarth)  transfor¬ 
mations  for  the  case  of  stagnation  flow 


T)  = 


X  Vw  i 1:3 

(Al) 

(du/dx)pex2  py 

7  1  pe  ^ 

(A2) 

in  which  the  velocity  at  the  boundary- layer  outer  edge  (inviscid  flow)  is 
written 
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Ue  =  X 


TJ  00 

B  —  x 


where  B  is  a  constant  and.  %  is  the  nose  radius.  A  similarity  parameter 
Ps  can  be  obtained  if  equation  (A2)  is  written  in  the  inverted  form 


R  R(du/dx)x2p( 


11  Pe 

fdT1 


Defining  Ps  as  the  coefficient  of  the  integral  and  substituting  equa¬ 
tion  (A3)  yields 

oor  < 

and  Ps  can  be  expressed  in  terms  of  Reynolds  number  based  on  wall  condi¬ 
tions,  Rew  =  PwRUoo/w  Using  the  equations  of  state  at  the  edge  and  wall 
with  py  =  pe  gives 

Pw  Pw 

Pe  ■  and  w  ‘  w 


Then  from  equation  (A5) 


_  l  l  zeTe 
n/2B  \/R eu  ZWTW 


Alternatively,  equation  (A5)  can  be  written  in  terms  of  the  stream  conditions. 
The  density  at  the  wall  is 

P  g(U  ) 


and  pe  p  .  Since  the  wall  temperature  and  U^,  are  assumed  to  be  constant, 


Py  -  Kl 


Writing  the  edge  density  as 


PooS 1  (  Uoo  )  PooToo 


6  "  ZeRTe 


g2(U  ) 
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and  substituting  this  and  p^-  from  equation  (A7)  in  equation  (A5),  we  have 

1 


Kn  Poo^W 


Z  T 

p  —  _ 6  e 

"s  n/2B  p  g2(U  )  /  U  RZttT 

00°^  '  00 '  v  00  W  c 


PS  =  g3(U00,T'oo) 


ZeTe 


M-W 


n/%  7  PmR 


(AS) 


If  the  effects  of  pressure  on  the  compressibility,  temperature ,  and  viscosity 
are  neglected,  the  similarity  parameter  from  equation  (A8)  is 


Rs  = 


F(Uoo^J 


(A9) 


As  in  the  case  of  two-dimensional  flow,  the  condition  for  similarity  of 
boundary- layer  solutions  in  the  axisymmetric  case  is  again  that  the  trans¬ 
formed  stream  function  at  the  wall  be  constant.  This  stream  function  is 


:W  "  -vw%  l  Pw^w  ^ 


du 


•1/2 


"  ”vwPw 


R 


XL 


>PWPWB, 


1/2 


(A10) 


To  show  that  this  condition  is  satisfied,  when  the  scale  varies,  we  obtain  an 
expression  for  fy  corresponding  to  the  ablation- product  flux.  In  the 
physical  system,  the  flux  is  v^py  =  m,  where  m  is  the  mass  rate  of  abla¬ 
tion,  in  mass  per  second  per  unit  area.  A  convenient  expression  for  m  (see, 
e.g.,  ref.  l6)  is 

m  =  4/Heff  (All) 


where  Kgfp  is  the  effective  heat  of  ablation,  assumed  to  be  a  function  of 
surface  material  and  stream  velocity,  and  q  is  the  cold-wall  convective 
heating  with  no  mass  injection.  The  assumption  that  the  heating  is  entirely 
convective  is  appropriate  for  the  velocities  considered  here,  but  for  some¬ 
what  higher  velocities,  effects  of  radiative  heating  would  also  have  to  be 
included.  Convective  heating  from  a  stagnation  region  may  be  obtained  from 
the  relation  given  in  reference  17, 


K2UVp  /R) 


1/2 


oo  oo 


(A12) 
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Equation  (All)  may  then  be  written 


“  =  \Pw=  ^/R  83  (Uj 


(A13) 


Expressing  stream  density  as 


P  = 
00 


and  vail  density 


Pfy  _  PooS(Uoo)  /^W 


and  substituting  equation  (A13)  in  equation  (AlO)  gives 

fw  =  -  / 


zwTw 


RRT, 


S3(Uoo)  J  PWU00B  PooHl(Uoo) 


=  -gJUoojTj 


(Aik) 


When  the  effect  of  pressure  on  compressibility  is  assumed  to  be  small  and 
the  viscosity  to  be  independent  of  pressure,  this  relation  shows  that  at 
constant  velocity,  the  transformed  stream  function  does  not  vary  with  body 
size  or  stream  density,  so  that  the  condition  for  boundary- layer  similarity 
is  met.1 

The  viscous  core  of  the  wake  is  considered  to  be  described  initially  as 
a  laminar  free- shear  layer  converging  to  a  minimum  diameter  at  a  constricted 
neck  located  about  one  body  diameter  downstream.  At  the  neck,  increased 
turbulence  is  assumed  to  disturb  the  laminar  temperature  and  velocity  profiles 
such  as  those  calculated  in  reference  19  so  that  temperature,  pressure,  and 
concentration  of  injected  gas  are  approximately  constant  across  the  diameter. 
Reference  4-  indicates  that  a  laminar  base  flow  with  turbulence  beginning  at 
the  neck  is  possible  throughout  a  large  range  of  stream  conditions.  On  the 
basis  of  experiments  described  in  reference  1,  and  the  analysis  in  refer¬ 
ence  19,  in  which  conditions  on  the  dividing  streamline  of  the  laminar  base 
flow  are  calculated,  the  velocity  of  the  flow  in  the  wake  core  is  assumed  to 
be  independent  of  scale. 

The  similarity  parameter  given  by  equation  (A9)  is  now  used  with  the 
assumptions  listed  above  to  estimate  effects  of  scale  on  radiation.  First, 
the  effect  of  varying  body  size  with  free- stream  properties  held  constant  is 
considered,  and  then  the  effect  of  maintaining  the  product  Rp^  constant 

^•This  result  is  also  shown  in  reference  l8  for  the  special  case  of  con¬ 
vective  heating  alone,  as  part  of  a  more  general  analysis  of  heating  at  higher 
speeds,  including  the  effects  of  radiative  heating  and  reradiation  from  heat 
shields. 
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(i.e.,  constant  Reynolds  number)  is  calculated.  The  first  approach  results  in 
a  very  simple  relation  for  the  scale  effect  on  radiation,  whereas  the  second 
approach  involves  a  more  extended  analysis  requiring  an  expression  for  radia¬ 
tive  intensity  as  a  function  of  gas  properties  and  concentration. 


Constant  Density 


Equation  (A4)  written  with  the  derived  expression  for  the  parameter  PQ 
(eq.  (A9))  becomes  s 

^  P. 


F(U  ) 
00 


dt] 

P 


(A15) 


When  the  stream  velocity  and  temperature  are  constant,  it  is  apparent  the 
boundary  layer  before  separating  has  a  thickness  inversely  proportional  to 
the  square  root  of  the  free-stream  Reynolds  number,  which  is  also  the 
classical  result  for  nonablating  surfaces. 

The  total  mass  flow  in  the  boundary  layer  (where  the  subscript  b  refers 
to  the  boundary- layer  edge  at  point  of  separation  from  the  body  surface)  is 
expressed  as 


where  5  =  y  at  the  point  where  up  «  U^.  If  the  similarity  parameter 

Ps  from  equation  (A9)  is  introduced  and  represents  the  constant  trans- 

formed  integral. 


%  KaR^PbUb^s^b  “  ^4Pbub  —  ^b  (Al6) 

n/^o 

Since  the  wake  core  is  considered  to  be  composed  of  the  boundary- layer  flow 
the  entrainment  of  air  due  to  mixing  in  the  separated  boundary  layer 
just  behind  the  body  is.  neglected),  the  mass  flow  in  the  wake  is  equated  to 
the  boundary- layer  flow 


^w  AwPw^w  “ 


(A17) 


where  the  subscript  w  refers  to  conditions  representing  the  viscous  wake  at 
the  station  at  which  radiation  effects  are  being  studied. 


The  mass  of  ablation  material  contained  in  the  wake  can  be  calculated 
from  the  ablation  rate  (eq.  (A13 ) )  and  the  effective  ablation  area  A^,  which 
is  proportional  to  R2.  The  rate  of  ablation  is  y 
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Mp  =  Apm 


=  K5R2  1-2.  ga(Uj 


(A18) 


Comparing  this  with  equation  (Al6) }  one  can  see  that  for  constant  veloc- 
ity  Uqo  and  constant  stream  density  p^,,  both  the  total  mass  flux  and  the 
mass  of  ablated  material  vary  with  R3/2.  Thus  the  relative  amounts  of  abla¬ 
tion  material  and  of  air  remain  constant  as  R  varies.  The  mass  of  these 
components  contained  in  a  length  Z  of  the  viscous  wake  can  be  computed  as 


iL 

Uw 


(A19) 


in  which  Uw  is  an  average  wake-flow  velocity  relative  to  the  body. 


With  equation  (A17) >  a  cross-section  area  and  radius  of  the  wake  can  be 
calculated: 


2  MB 

A»  "  ■*»  =  ZJi. 


(A20) 


w  w 


If  we  substitute  Mp  from  equation  (Al6)  and  solve  for  Rw 


*w  =  Ke| 


K 


3/2 


PbUb  _ _ 

Pw^w 


1/2 


(A21) 


It  is  assumed  that  the  pressures  pb  at  the  body  base  and  pw  in  the 
wake  are  constant  fractions  of  the  stream  pressure.  Then 


Pb  = 


ZbRTb 


ZbTb 


(A22) 


and 


w 


Pv  Vw 


ZWTW 


(A23) 


In  these  expressions  Up,  Tp,  and  Zp  are  determined  by  the  inviscid 
flow  and  do  not  vary  with  scale.  As  mentioned  earlier,  the  wake  velocity  is 
relatively  independent  of  scale.  Equation  (A12)  then  is 


,3/2 


1/2 

Rw  =  MW*""  1  (A24) 

The  wake  temperature  Tw  varies  to  some  extent  as  a  result  of  variations  in 
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the  flow  deflection  at  the  base  of  the  body  and  at  the  re compression  shock 
wake;  however,  estimates  indicate  that  the  variation  behind  the  trailing 
shock  is  not  large.  If  it  is  assumed  that  this  variation  and  the  variation 
in  Zw  can  be  neglected,  equation  (A24)  reduces  to 

(A25) 

This  analysis  indicates  that  a  variation  in  scale  at  constant- stream 
density  and  velocity  produces  no  variation  in  the  absolute  densities  of  the 
two  constituent  gases,  but  a  variation  in  wake  radius  given  by  equation  (A25) . 
On  the  supposition  that  the  radiative  intensity  is  determined  by  the  concen¬ 
tration  of  the  injected  gas  and  the  local  properties  of  the  mixture,  the 
radiation  from  a  fixed  volume  (watts/cm3)  is  independent  of  scale.  A  possible 
restriction  on  the  use  of  this  result  arises  because  of  the  assumption  that 
the  wake  is  optically  thin.  Among  the  materials  studied,  those  producing 
bright  wakes  may  be  well  represented  by  this  assumption  only  at  the  lower 
densities  and  smaller  scale  factors. 


3  =  KsR-l/4 


Constant  Reynolds  Number 


The  effect  on  the  peak  radiative  intensity  of  varying  scale  while 
maintaining  constant  Reynolds  number  (based  on  free- stream  conditions)  is 
investigated  in  this  section.  Constant  Reynolds  number  here  is  meant  to 
express  p^R  =  constant.  Equation  (A2l)  can  be  written 


R. 


w 


_  M)ub 


R 


PvUv 


(Bej'l/4(lb)1,a 


(A26) 


As  before,  using  the  assumption  that  the  velocities  and  Uv  are  constant 
fractions  of  U^,  and  inserting  equations  (A22)  and  (A23)  in  equation  (A26) 
gives 


1^, 

N/ZbTb 


(A27) 


In  the  approximate  analysis  presented  here,  it  is  consistent  to  assume  again 
that  variations  in  temperature  and  compressibility  ratios  can  be  neglected. 
Then 

h  =  K9(Reoo)"1/4  (A28) 


The  variation  of  wake  radius  given  by  equation  (A28)  has  been  compared  with 
that  obtained  from  the  tests.  The  experimental  results  shown  in  figure  l6 
were  obtained  by  measuring  the  width  of  the  luminous  wake  recorded  with  image 
converter  and  Kerr-cell  camera  photographs.  The  one-fourth  power  variation 
is  in  agreement  with  the  measurements. 
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For  scale  variations  with  the  product  p^R  constant,  the  flow  geometry 
and  transition  characteristics  are  unchanged.  Also,  from  the  boundary- layer 
similarity  conclusions,  the  ablation  gas  concentration  does  not  vary.  This 
is  shown  by  the  ratio  %/%  where  ftp  is  given  by  equation  (Al6)  and  Mp 
is  obtained  from  equation  (Al8) .  The  mass  ratio  then  is 


— =-  =  K- 


p  g  (U  ) 

oo  2  00 

”p,  uTiT" 

b  b 


(A29) 


If  the  effects  of  scale  on  the  temperature  and  compressibility  are  again 
assumed  to  be  negligible ,  this  ratio  is  a  constant. 

Radiant  intensity  equation.-  With  p^R  constant,  changes  in  scale  lead 
to  inverse  changes  in  density.  Thus,  an  expression  is  needed  to  relate 
observed  radiant  intensity  to  that  at  conditions  to  which  the  scaling  rela¬ 
tions  apply.  The  following  relation  was  used  in  reference  8  to  express 
radiative  output  from  ablation  products  in  a  boundary  layer.  If  I  is  the 
power  per  unit  volume, 


i/pp  =  rp^T|  ,  watts/g 


(A30) 


where  pp  is  the  density  of  injected  ablation  gas,  p^  is  the  total  gas 
density,  F  is  a  constant, and  £,  and  |  are  exponents  to  be  evaluated.  Apply¬ 
ing  this  formula  to  radiation  per  unit  length  of  wake  gives 


AWI  =  EAwppp^4  >  watts /cm 


(A31) 


Noting  that  pp  =  mAp/uwAw  where  Ap  is  the  effective  ablation  area,  equa¬ 
tion  (A13)  furnishes  a  relation  for  the  ablation  gas  density 


Pp  = 


3/2  1/2  ,TT  « 

R  p  g2  (U  ) 

*00  x  oo7 


UA- 


(A32) 


Combining  equations  (A17)  and  (Al6)  gives  the  total  density  in  the  wake,  Pw, 


K4  _  tt  o3/2„  -1/2t 

Pw  _  A  tt  PbUbR  Pco  Ib 


AWUW 


(A33) 


When  the  proportionalities  U’0  cc  and  oc  p^/Z^Th  are  inserted,  equa¬ 
tion  (A31)  is 


A  J  vp3/2  1/2  g2(U»YR  IbPbo  Uoot  rpi 

AyX  -  KR  v 


(A3*0 
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Making  use  of  the  previously  listed,  assumptions  (Z^,  T^,  and  Uw 
constant)  in  equation  (A34),  we  have 


AWI  =  K12R  p^T| 


(A35) 


If  the  following  is  substituted  for  wake  area  (from  eq.  (A28)) 


A 


w 


equation  (A35)  gives  the  following  radiant  power  per  unit  volume. 


I  = 


KisrEt1 

R 


(A36) 


(A37) 


Thus  except  for  the  nonzero  exponent  on  pro  and  the  variation  of  temperature 
with  density ,  the  optical  thickness  is  preserved  if  Reynolds  number  is  held 
constant j  that  is ,  the  emissive  power  is  inversely  proportional  to  the  scale 
factor.  If  the  test  velocity  is  adjusted  to  hold  temperature  constant  and 
minor  adjustments  of  density  are  made  to  preserve  Reynolds  number,  better 
simulation  of  large-scale  bodies  might  be  achieved  in  the  ballistic  range. 

If  we  insert  the  observed  variation  with  density  of  radiant  power  per 
body  diameter  (let  |  take  on  the  value  4  (ref.  8))  and  include  the  varia¬ 
tion  of  temperature  with  density  calculated  for  equilibrium  air  at  constant 
enthalpy,  the  following  approximate  expression  is  obtained: 


I  =  Kl4R'3/2 


or 


W  =  2*R^RI  ex 

v/^e 


(A38) 

(A39) 


w  =  k15r3/2 


(A40) 


This  furnishes  a  relation  for  scaling  the  radiant  power  from  a  length  of  wake 
equal  to  the  body  diameter  at  constant  Reynolds  number. 
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Figure  1.-  Schematic  drawing  of  pilot  hypersonic  free-f light  facility. 


NOTE:  Dimensions  in  cm 


Figure  3.-  Radiometer  and  slit  arrangement. 
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Oscillograms  shoving  radiometer  output  signals  at  a  series  of  wavelength; 
Uqo  =  6.6  km/sec,  P^/Pq  =  0.07,  polycarbonate. 


b)  Image  converter  photograph  of  luminous  flow;  polycarbonate  models 
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( a )  Polyc  arb  onat  e . 

Figure  7* -  Peak  wake  spectra  at  various  stream  densities 


Wavelength,  X,  microns 

(b )  Polyethylene • 
Figure  7*“  Concluded. 


□  □ 


*e  8.-  Measured  radiative  power  (in  spectral  region  from  0.2  to  1.0  p) 
wake  at  location  of  peak  radiance  as  a  function  of  stream  density. 


Radiant  power,  W,  watts/diam 


Velocity,  Uoo,  km/sec 

(b )  Polyethylene . 
Figure  Q* ”  Concluded 


various  longitudinal  locations  behind  model. 


Figure  10.-  Continued 


Pco/Po  ~  0-018 
Uoo  =  6.9  km/sec 


(c)  G.  E.  124. 


Figure  10.-  Concluded. 


L^eouy  L/U“  i 


Distance,  x/D 

igure  12.-  Longitudinal  distribution  of  wake  luminosity  in  the  wavelengtl 
region  0.45  to  0.59  micron  for  polyethylene  models. 


A-30236 

Figure  l4 .  -  Successive  photographs  of  luminous  flea-/  surrounding  and  in  wake 
of  polyethylene  model;  Re^  =  68,400;  Uro  =  7.8  km/sec;  pbo/p0  =  0.02. 
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Ambient  density,  atm 


Figure  15.-  Ratio  of  total  wake  radiative  power  to  total  power  from  equilibrium 
air  in  the  model  gas-cap  polycarbonate;  D  =  0.71  cm. 


Figure  l6,-  Minimum  diameters  of  luminous  wakes 
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